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ABSTRACT 

Context. The apparent position of the "core" in a parsec-scale radio jet (a compact, bright emitting region at the narrow end of the jet) 
depends on the observing frequency, owing to synchrotron self-absorption and external absorption. While providing a tool probing 
physical conditions in the vicinity of the core, this dependency poses problems for astrometric studies using compact radio sources. 
Aims. We investigated the frequency-dependent shift in the positions of the cores (core shift) observed with very long baseline 
interferometry (VLB!) in parsec-scale jets. We discuss related physics, as well as its effect on radio astrometry and the connection 
between radio and optical positions of astrometric reference objects. 

Mettiods. We searched for the core shift in a sample of 277 radio sources imaged at 2.3 GHz (13 cm) and 8.6 GHz (4 cm) frequency 
bands using VLBI observations made in 2002 and 2003. The core shift was measured by referencing the core position to optically 
thin jet features whose positions are not expected to change with frequency. 

Results. We present here results for 29 selected active galactic nuclei (AGN) with bright distinct VLBI jet features that can be used in 
differential measurements and that allow robust measurements of the shift to be made. In these AGN, the magnitude of the measured 
core shift between 2.3 and 8.6 GHz reaches 1.4 mas, with a median value for the sample of 0.44 mas. Nuclear flares result in temporal 
variability of the shift. 

Conclusions. An average shift between the radio (4 cm) and optical (6000 A) bands is estimated to be approximately 0. 1 mas, and 
it should be taken into account in order to provide the required accuracy of the radio-optical reference frame connection. This can be 
accomplished with multi-frequency VLBI measurements yielding estimates of the core shift in the sources used for the radio reference 
frame and radio-optical position alignment. 

Key words, galaxies: active — galaxies: jets — radio continuum: galaxies — astrometry — reference systems 

1 . Introduction frequencies can be used to determine critical physical and geo- 
metrical parameters of the relativistic jet origin. 
Extragalactic relativistic jets are formed in the immediate vicin- The core shift induced by synchrotron self-absorption ex- 
ity of the central black holes in galaxies, at distances on the order tends over the entire spectral range and has an immediate ef- 
C3 of 100 gravitational radii, and they become v isible in the radio at feet on astrometric measurements made in the radio and op- 
distan ces of about 1000 gravitational radii (iLobanov & ZensusI tical domains. One particular issue likely to be affected by 
|2007|). This apparent origin of the radio jets is commonly called the core shift is the connection between the radio and opti- 
the "core". In radio images of extragalactic jets, the core is lo- cal positions of distant quasars. At present, the International 
cated in the region with an optical depth T, * 1 . This causes the Celestial Reference System (ICRS) is realized prii narily by 
absolute position of the core, r^o^, to vary with the observing fre- the International Cele stial Reference Frame (ICRF, iMa et alj 
quency, v, since the opt ical depth profile alon g the jet depends on |1998|: iFev et al.l l2004l) based on very long basehne interfer- 
v: rcoie oc v^'/*' (Blan dford & Konigl| [T979l). See an illustration ometry (VLBI) measurements at 2.3 GHz (13 cm band) and 
ofthis eff^ect in Fig.[T] Variations in the optical depth along the jet 8.6 GHz (3.6 cm band). The present extension of the ICRF 
can result from synchrotron self-absorption (Konigl 1981), pres- to visible Ught is based on the Hipparcos Catalogue, with 
sure and density gradients in the jet and free-free absorption in rms uncertainties estimated to be 0.25 mas yr^^ in each com- 
the ambient m edium most likel y associated with the broad-line ponent of the spin vector of the frame and 0.6 mas in the 
region (BLR) (lLobanovl lI998bl). If the core is self-absorbed and compone nts of the orientation v ector at the catalogue epoch, 
in equipartition, the power index k, = 1 dBIandford & Konigll J1991.25 jKovalevsky et alj|1997 '). The link between the VLBI 
[1979). Density and pressure gradients in the je t and external ab - an d Hipparcos referenc e frames and its accuracy is discussed 
sorption can lead to deviations in k, from unity (lLobanovlI998bh . by iLestrade et all (Il995h . In the near future, new optical astro- 
Changes in the core position measured between three or more metric reference frames will be established by the GAIA as- 
trometry mission (Lindegren& Ferryman 1996) and the Space 
Interferometry Mission (^/M. lShaoll99 8). Both missions will ex- 
Send offprint requests to: Y. Y. Kovalev plore largely the 3000 A-10000 A range coiTesponding to (0.3- 
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Fig. 1. A scheme illustrating the frequency-dependent position shift of the VLBI core. Adopted from lLobanovl (1 1 996h . 
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Fig. 2. Example of global VLBI results. 2.3 and 8.6 GHz naturally weighted images of the quasar J1642+6856, epoch 6 March, 
2002. The beam, with typical full width at half maximum of 2.6 mas at 2.3 GHz and 0.7 mas at 8.6 GHz, is shown in the lower left 
hand corner of each map. The contours are plotted in successive powers of 2 times the lowest contour. The peak intensity and the 
lowest contour for the 2.3 GHz and 8.6 GHz images are 367 and 1.2 mJybeam 752 and 0.8 mJy beam respectively. The core 
('C') position is measured relative to the optically thin jet feature 'A'. The spectral index of the core is ff = 0.6, while that of the jet 
feature a - -0.6 {S oc v"). 



1) X 10'' GHz. GAIA will measure optical positions of about VLB AQ CaUbrator Survey ([Beaslev et al.ll2002l: iFomalont et"al 



500,000 distant quasars with < 0.1 mas accuracy, similar to 2003; Petrov et al. 2005, 2006", Kovalev et al. 2007; Petro v et al. 



the accuracy of radio source positions from the ICRF and the i2008.) . SIM is expected to deliver high precision astrometric po- 
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sitions of bright quasars at ^ lOjuas accuracy and a target ac- 
curacy of a: 20 //as is envisaged for the radio-optical alignment 
jUnwin 2005). For both of these missions, an important prob- 
lem will be matching the optical catalogues to the radio as- 
trometry catalogues based on precise positions of compact ex- 
tragalactic objects and determini ng mutual rotations, distortions 
and zonal systematic errors (e.g. Fey et al. 2001; Souchay et al.' 
[2006; Lambert et al. 2006; Frey et al. 2006; Bourda et al. 2007). 
This match relies on an assumption that the dominating compo- 
nent of emission in both the radio and optical bands is physically 
the same region (see sect. l3.4l for detailed discussion). 

The core shift is expected to introduce systematic offsets be- 
tween the radio and optical positions of reference sources, af- 
fecting strongly the accuracy of the radio-optical matching of 
the astrometric catalogues. The magnitude of the core shift can 
exceed the inflated errors of the radio and optical positional mea- 
surements by a large factor This makes it necessary to perform 
systematic studies of the core shift in the astrometric samples in 
order to understand and remove the contribution of the core shift 
to the errors of the radio-optical position alignment. 

Measurements of the core shift have been done so far 
only in a small number of objects (e.g ., [Marcaide et al. 1994: 
Lara et al."1994'; 'Porcas & Rioia'4997^ Cobanov 1996, 1998b: 
Pai-asi et al. 2000; Ros & Lobanov 2001; Bietenholz et al. 2004: 
Kadler et alj i2004). In this paper, we present for the first time 



results for 29 compact extragalactic radio sources used in VLBI 
astrometric studies. In Sect. |2] we describe global VLBI data, 
the source sample, and the method of measurement used in this 
paper, and present results of the core shift measurements. In 
Sect. [3] we discuss astrophysical applications of the frequency- 
dependent core shift (Sect. 13.11 ). variability of the shift result- 
ing from nuclear flares (Sect. 13.21 ), effect of the shift on multi- 
frequency VLBI studies (Sect. [373T i. We investigate the influence 
of this effect on the radio-optical reference frame matching and 
suggest a method to compensate for it (Sect. |3.4T l. We summarize 
our results in Sect. |4] 



2. Core shift measurements 

We have imaged and analyzed 277 sources from geodetic 
observations made in 2002 and 2003 — ten 24hr-long ex- 
periments on 16 January 2002, 6 March 2002, 8 May 2002, 
24 July 2002, 25 September 2002, 11 December 2002, 
12 March 2003, 7 May 2003, 18 June 2003. Geodetic RDV ses- 
sions feature simultaneous observations at 2.3 GHz and 8.6 GHz 
(S and X bands) with a global VLBI network at right circular 
polarization. This includes for every session the VLBA and up 
to nine other radio telescopes from the following list: Algonquin 
Pai-k (46 m), Gilcreek (26 m), HartRAO (26 m), Kokee (20 m), 
Matera (20 m), Medicina (32 m), Noto (32 m), Ny Alesund 
(20 m), Onsala (20 m), TIGO (6 m), Tsukuba (32 m), Westford 
(18 m), Wettzell (20 m). The data processing technique and 
imaging results are described by Pushkarev & Kovalev (2008). 

This long-term RDV program is one of the best choices for 
a large project to measure two-frequency core shifts on the basis 
of open archival raw VLBI data for several reasons: (i) it is op- 
timized to have a good (M,v)-coverage, (ii) it has the maximum 
possible resolution for ground-based VLBI at these frequencies, 
(iii) the frequency ratio between the simultaneously observed 
bands is high (3.7), and (iv) the core shift per unit of frequency 



between 2.3 and 8.6 GHz is larger th an that at higher frequencies 
because of opacity effects (see, e.g..'Loban ovll998bl) . 

We have measured the frequency-dependent core shift 
by model-fitting the source s tructure with two-dimensional 
Gaussian components ( IPearsonllT999i) and referencing the posi- 
tion of the core component to one or more jet features, assuming 
the latter to be optically thin and having frequency-independent 
peak positions (Fig.[TJ. The method is illustrated inFig.|2] where 
astrometric S/X images of 1642-1-690 are shown, with positions 
marked for the VLBI core ('C') and jet feature ('A'). 

D. C. Homan and Y. Y. Kovalev (in prep.) have made tests 
comparing core shifts measured in the quasars 1655H-077 and 
2201-1-315 with relative astrometry (phase referencing to a cali- 
brator source) to those obtained as here from self-calibrated im- 
ages by referencing the core to optically thin jet features. The 
core shifts obtained using these two methods agreed within the 
errors. This result indicates that self-calibrated images provide 
sufficient accuracy for measuring the core shift and can be used 
for a systematic study of this effect. However, the VLBI phase 
referencing method is required for jets whose structure does not 
allow the self-referencing technique to be applied. 

Other approaches have also been tried. iRioia et akl (|2005|) 
have proposed a "source/frequency phase referencing" method 
for measuring the core shift, which could be particularly useful 
while working at higher radio frequencies. I Walker et al.l (l2000h 
performed a 2-D cross-correlation analysis in order to align 
VLBA images of 3C 84 measured at different frequencies with 
a formal a ccuracy less than 5 uas, s ee also discussion of this 
method by ICroke & Gabuzdal ( 20081) . However, the real error 
margins of the cross-correlation method are certainly substan- 
tially larger, as this method requires the assumption that there 
are no large variations in observed spectral index along the jet 
(if the core region is blended) which is not the case for many 
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Fig. 3. Histogram of the derived core shift values for 29 sources. 
One average core shift value per source is used. The median 
value for the distribution is equal to 0.44 mas. 
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Table 1. Results of core shifts measurements for 29 bright extragalactic radio sources 



lAU name 


J2000.0 
R.A. 




J2000 
Decl. 




redshift 


Observing 
epoch 


Core shift" (mas 
between 2.3 and 8.6 


0003-066 


00 


06 


13 


893 


-06 


23 


35 


335 


0.347 


2003-05-07 


+0.375 


+ 0.019 






















2003-09-17 


+0.246 


± 0.028 


0118-272 


01 


20 


31 


663 


-27 


01 


24 


652 


>0.557 


2003-03-12 


+0.791 


± 0.075 


0148+274 


01 


51 


27 


146 


+27 


44 


41 


794 


1.26 


2002-09-25 


+ 1.229 


+ 0.018 


0202+149 


02 


04 


50 


414 


+ 15 


14 


11 


044 


0.405 


2002-03-06 


+0.058 


± 0.043 


0239+108 


02 


42 


29 


171 


+ 11 


01 


00 


728 


2.680 


2003-06-18 


+0.868 


+ 0.032 


0342+147 


03 


45 


06 


417 


+ 14 


53 


49 


558 


1.556 


2002-01-16 


+0.274 


+ 0.034 


0425+048 


04 


27 


47 


571 


+04 


57 


08 


326 


0.517 


2002-09-25 


+0.385 


+ 0.046 


0430+052 


04 


33 


11 


096 


+05 


21 


15 


619 


0.033 


2003-06-18 


+ 1.196 


+ 0.018 






















2003-09-17 


+0.956 


±0.016 


0507+179 


05 


10 


02 


369 


+ 18 


00 


41 


582 


0.416 


2002-07-24 


+0.383 


±0.012 


0607-157 


06 


09 


40 


950 


-15 


42 


40 


673 


0.324 


2002-01-16 


+0.140 


± 0.038 






















2002-01-16 


+0.149 


± 0.065 


0610+260 


06 


13 


50 


139 


+26 


04 


36 


720 


0.580 


2002-12-11 


+ 1.362 


± 0.030 






















2002-12-11 


+ 1.370 


± 0.030 


0736+017 


07 


39 


18 


034 


+01 


37 


04 


618 


0.191 


2003-05-07 


+0.186 


± 0.083 


0839+187 


08 


42 


05 


094 


+ 18 


35 


40 


991 


1.272 


2003-03-12 


+0.911 


±0.017 


0850+581'' 


08 


54 


41 


996 


+57 


57 


29 


939 


1.322 


1997-01-10 


+ 1.349 


±0.031 


0952+179 


09 


54 


56 


824 


+ 17 


43 


31 


222 


1.478 


2002-01-16 


+0.725 


± 0.039 


1004+141 


10 


07 


41 


498 


+ 13 


56 


29 


601 


2.707 


2002-01-16 


+0.399 


± 0.054 






















2003-05-07 


+0.490 


± 0.034 


1049+215 


10 


51 


48 


789 


+21 


19 


52 


314 


1.300 


2002-01-16 


+0.997 


± 0.068 


1 147+245 


11 


50 


19 


212 


+24 


17 


53 


835 


0.200 


2003-05-07 


+0.321 


±0.017 


1219+285 


12 


21 


31 


691 


+28 


13 


58 


500 


0.102 


2002-03-06 


+0.239 


± 0.008 


1642+690 


16 


42 


07 


849 


+68 


56 


39 


756 


0.751 


2002-01-16 


+0.475 


± 0.024 






















2002-03-06 


+0.437 


± 0.025 






















2003-05-07 


+0.350 


± 0.039 


1656+053 


16 


58 


33 


447 


+05 


15 


16 


444 


0.879 


2002-09-25 


-0.010 


± 0.048 


1655+077 


16 


58 


09 


Oil 


+07 


41 


27 


541 


0.621 


2002-07-24 


+0.672 


± 0.049 


1803+784 


18 


00 


45 


684 


+78 


28 


04 


018 


0.680 


2002-05-08 


+0.722 


± 0.022 






















2002-05-08 


+0.645 


± 0.053 


1830+285 


18 


32 


50 


186 


+28 


33 


35 


955 


0.594 


2002-03-06 


+0.772 


±0.012 


1845+797 


18 


42 


08 


990 


+79 


46 


17 


128 


0.057 


2002-01-16 


+0.551 


± 0.049 






















2002-01-16 


+0.403 


± 0.083 


1926+087 


19 


28 


40 


855 


+08 


48 


48 


413 




2002-03-06 


-0.088 


± 0.024 


2021+317 


20 


23 


19 


017 


+31 


53 


02 


306 




2002-07-24 


+0.158 


±0.105 


2143-156 


21 


46 


22 


979 


-15 


25 


43 


886 


0.698 


2002-05-08 


+0.121 


±0.018 


2155-152 


21 


58 


06 


282 


-15 


01 


09 


328 


0.672 


2003-09-17 


+0.369 


±0.018 



Column designation: Col. 1 - lAU source name (B1950), Col. 2 - Right Ascension (J2000), Col. 3 - Declination (J2000), Col. 4 - redshift as 
given by Veron-Cettv & Veron (2006) with the following exceptions: the redshifts for 0239+ 10 8, 0342+147, 1 147+245, 2320+506 are taken from 
ISow ards-Emmerd et al. (2005); the redshift for 0425+048 is taken from Afanasiev et alj 12003"), Col. 5 — epoch at which the shift was measured 
(YYYY-MM-DD), Col. 6 - shift of the VLBI core position (mas) between 2.3 and 8.6 GHz and its Icr error. If the shift is measured relative to two 
different jet features, two independent values for the same epoch are given. 

" Positive values of the shift mean that the distance between the core and the jet feature measured at 8.6 GHz is greater than the one measured at 
2.3 GHz. 

* The measurement for this source is based on our imaging results of the NRAO archival data from the VLBA only S/X project BF025 (see 
description of the BF025 program in lFev & Charlotl[2000) . 



extragalactic jets (see, e.g., the reconstructed spectral index gra- 
dient along the jet for 3C 120 in Fig.|7]l. 

We have measured core shifts between 2.3 and 8.6 GHz in 
29 AGN (see Fig.[3]and TablelTJ, with the resulting values of the 
shift ranging between -0.1 and 1.4 mas and the median value 
for the sample of 0.44 mas. Errors of the core shift measurements 
presented in Table[T]are estimated fro m the uncer tainties of com- 
ponent positions calculated following lFomaloni ([T999i) . The ro- 
bustness of the measurements is supported by the following ob- 
servation. In objects with several bright jet components, similar 
core shift values within errors are obtained from referencing the 
core position to different jet features. The core shift measure- 
ments for these objects are indicated by separate symbols in the 



top panel of Fig. |4] In objects observed more than once, com- 
parable core shift values are obtained from referencing the core 
position to the same jet feature at different observing epochs (see 
sect. 13.21 for details on a possible core shift variability). In both 
cases, no systematic dependencies are observed (Fig.|4|, imply- 
ing that measured core shift values do not depend on the compo- 
nent location and features cross-identified between the two fre- 
quency bands S/X have similar position angles (P.A.f^^p). In or- 
der to further ensure the fidelity of the measurements made, we 
have also performed an independent check by analyzing spectral 
index images with and without core shift applied and dropping 
all objects in which application of the core shift resulted in un- 
physical values of the spectral index (see details below). 
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Fig. 4. Dependencies between several derived parameters to 
check for possible systematics presented in results. 
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Fig. 5. Distribution of blending shifts of the core position. The 
distribution is calculated for all the sources and epochs in Table[T| 
at 8.6 GHz by downgrading the resolution to the one at 2.3 GHz. 
Positive values correspond to inward shifts of the core. 



The difference in resolution and corresponding angular size 
of the restoring VLBI beam at 2.3 GHz and 8.6 GHz could lead 
to a systematical apparent (not real) shift of the core position 
due to resolution-dependent blending of the core and optically 
thin parts of the jet. To estimate the magnitude of this effect, we 
have downgraded the resolution of 8.6 GHz data to the one of the 
2.3 GHz band by setting the maximum (M,v)-radius at 8.6 GHz 
to about 75 MA. Then the 8.6 GHz model was re-fitted for every 
source and epoch from Table [1] by varying all free parameters. 
These models were compared to the original 8.6 GHz ones, and 
the magnitude of the core shift due to blending was calculated. 
Distribution of these blending shifts of the core position along 
the jet axis at 8.6 GHz is presented in Figure |5] The median 
value for this blending effect is ^ 0.006 mas which is much less 
than the typical error of the core-shift measurements presented 
in this paper (Table [T]). 

For about 90 percent of the 277 objects imaged, no reliable 
estimates of core shifts have been obtained. All of these objects 
satisfy one or more of the following five conditions: 



1) The core and jet features could not be extracted reliably at 
one or both frequency bands because of components blend- 
ing, or a smooth or weak jet. 

2) The same structures at both bands could not be identified 
reliably. 

3) Modeling by two-dimensional Gaussians is not adequate due 
to the complexity of the observed brightness distribution. 
Decision on every of the items 1) to 3) is made on the ba- 
sis of our analysis of two-frequency model parameters and 
their errors, values of fits of models to calibrated visibil- 
ity datasets, as well as comparison of residu al images after a 
model is subtracted with the imaging noise ( iFomalontll 19991: 
lPearson|[T999i) . 
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4) Jet components may be partially opaque (spectral index a > 
0, 5 oc v"^) and as a result of this their positions change with 
frequency. 

5) Application of derived core shift value results in unphysical 
values of the spectral index. This implies large optical depth 
observed in extended jet regions (a > 0), and/or optically 
thin core region (a < 0), and/or steep and irregular a gradi- 
ents across the jet. 

The goal of this study is to provide and analyze most reliable 
core shift detections and therefore such stringent selection cri- 
teria are applied to reduce the chance of an erroneous measure- 
ment. Based on these considerations, we deem acceptable the 
relatively low percentage (about 10 percent of the sample) of 
positive measurements. A designated study with significantly 
higher rate of successful core shift detections would require 
dedicated, time consuming, multi-band VLBI experiments with 
phase referencing. We estimate that about 6 hours of VLBA time 
per source will be needed in this case, assuming the current sus- 
tainable recording data rate of 128 Mbps. The time requirements 
would be shorten significantly after the VLBA recording system 
is upgraded for operations at a 4 Gbps recording rate. 

3. Physics and astrometry with the core shift 

The effect of frequency-dependence of the core position has an 
immediate connection to several physical and astrometric studies 
using compact extragalactic radio sources. Systematic observa- 
tions of the core shift in a sample of compact jets can be used for 
probing the conditions in the compact jets and nuclear regions 
in AGN and understanding the effect the shift may have on the 
alignment of the radio and optical reference frames. 

3. 1 . Deriving piiysical properties of tiie opaque jet base and 
absorbing material 

If the core shifts and the power index are measured from 
VLBI observations at more than two frequencies, the magnetic 
field strength and distributio n can be recon structed in the ultra- 
compact regions of the jets dLobanovll 1998 b). The offset of the 
observed core positions from the true base of the jet as well as 
the dist ance from the n ucleus to the jet origin can also be de- 
rived dLobanovll 1998bl) . Estimates of the total (kinetic + mag- 
netic field) power, synchrotron luminosity and the maximum 
brightness temperature, Tb^max, in the jets can be made. The ra- 
tio of particle energy to magnetic field energy can also be esti- 
mate d, from the de rived rb,max- This information enables testing 
of the Konigl ( 198 1 ) model and several of its later modifications 
(e.g., Hutter & Mufson 1986; Bloom & Marscher 1996). The es- 
timated distance from the nucleus to th e jet origin constrains the 
self-si milar jet model (iMarscheill 1 9951) and the particle-cascade 
model (iBlandford & Levinsonl ll995). This approach can also 
be applied to determine the matter content in parsec-scale jets 
jHirotani et al. 2000 ; Hiro tani 20051) . 

If the core shifts are measured at four or more frequencies, 
the absolute geometry of the jet can be determined, giving the 
absolute offset, rcoie [pc], of the core from the central engine. 
Combined with the value for magnetic field, Score [G], this gives 
an estimate for the mass of the central black hole Mbh[M0] w 

7 X lO^Beore'-core (lLobanovll998bl) . 

Density and pressure gradients in the jet and external absorp- 
tion can lead to deviations in ki- from unity. If the changes in ki- 
with frequency are measured with sufficient precision, they can 
be used to estimate the size, particle density and temperature of 



the absorbing material surrounding the jets dLobanovll 1998bl) . 
These estimates can be compared with the black hole masses 
and BLR sizes obtained from reverberation mapping and appli- 
cations of the relation established between the black hole mass 
and various ki nematic properties of central regions in AGN (c.f., 
iPeterson & Be ntz 2006). 

3.2. Nuclear flares and core shift 

Opacity in the compact jet can be changed substantially if the 
source undergoes a nuclear flare. This can affect the observed 
position of the core at a given frequency and modify the core 
shift between different frequencies. The magnitude of this effect 
can be estimated using the following argument. Let us assume 
that, during a flare, the compact jet does not change its orienta- 
tion and opening angle. Then, following iLobanov, (.1998b.) . the 
magnitude of the core shift is 



2/3 ^2/3g-2/3p-4/3 



(1) 



Fi, 6i are the 



where 5 core is the flux density of the core and y6j, ij, 
speed, Lorentz factor and Doppler factor of the jet, respectively. 
The flux density of the core is 5 core °^ "Mj^coi"' where Score 



IS 



the magnetic field in the core region, A^o is the particle density of 
the jet plasma and a is the spectral index. It follows immediately 
that 



(2) 



It has been shown that nuclear flares are likely to be pro- 
duced by particle density variations in the compact jet, while 
the Doppler factor and magnetic field change only weakly 
(iLobanov & Zensusi 1 19991) . Thus, it can be expected that, in 
the simplest case Ar oc NfJ^ oc 5core- This proportionality 



m 

(U 

o 




8.6 GHz core flux ratio 

Fig. 6. Core shift ratio versus 8.6 GHz core flux density ratio. 
The ratios are calculated for the four sources in the sample which 
have measurements at more than one observing epoch (Table[T|i. 
The core shift ratio is defined to be greater than unity, i.e., for 
each pair of measurements, the larger value of the shift is in the 
numerator of the ratio. The line is a least square fit to the points, 
yielding a slope of 0.96 + 0.18. 
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Fig. 7. Influence of the core shift on positional alignment of multi-frequency VLBI images, illustrated by 2-8 GHz spectral index, 
a (S oc v"), images for the jet in 3C 120. Global VLBI observations at 2.3 and 8.6 GHz made on 18 June 2003, were used. In each 
panel, contours represent the total intensity distribution at 8.6 GHz. The restoring beam used for convolution of the total intensity 
images at both frequencies is shown in the lower left hand corner of each map. The spectral index image on the top left is made by 
aligning the 2.3 and 8.6 GHz images at their respective phase centers. The top right image is aligned at model-fitted core positions. 
The bottom image is produced by aligning optically thin jet components (the arrow marks the resulting core shift between the 2.3 
and 8.6 GHz images). 



can be measured from comparison of core shifts measured in 
the same object at different epochs, yielding a power index 
e = log(Ar,/Ar2)/log(5,/52). 

In our data, core shifts can be compared at multiple epochs 
only in four objects, at five instances altogether (see Table [1]). 
The ratios of the core shifts measured at different epochs in the 
same objects are plotted in Fig.|6]against the respective ratios of 
the core flux density at 8.6 GHz. Ideally, the ratios of core shifts 
should approach unity as the flux density ratio approaches unity. 
This is not the case in Fig. |6] implying that flares in compact 
jets involve changes not only in the density of the plasma but 



also in its physical properties. This can be further examined by 
analyzing the slope of the relation plotted in Fig.|6] 

The measurements plotted in Fig.|6]indicate a positive trend, 
with larger core shifts for higher core flux density. A linear fit 
shown in Fig.|6]corresponds to e = 0.96 + 0.18, which is larger 
than the value of 2/3 expected for nuclear flares resulting entirely 
from plasma density variations in the flow. This may imply that 
flaring jets contain a hotter plasma or stronger magnetic field. 
However, the small number of data points in the present plot re- 
sults in large errors and precludes making any conclusions about 
the physical meaning of the fit, and we present it here only as an 
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example of the effect nuclear flares may have on the core shift. 
This effect certainly deserves a more comprehensive study. 



3.3. Effect of the core shift on multi-frequency VLBI studies 

In VLBI observations, application of the self-calibration 
technique based on closure phases (e.g., IjennisonI Il958t 
iReadhead & WiUdnsonI flOTSh results in a loss of information 
about absolute positions of target objects. In this situation, 
multi-frequency VLBI images are often aligned at the appar- 
ent position of the core at different frequencies (i.e., assum- 
ing implicitly a zero core shift). In the case of a non-zero core 
shift, this approach may undermine the results of any study 
based on multi-frequency VLBI data (e.g. spectral imaging and 
Faraday rotation measurements). It has been shown that ac- 
cou nting for the core shift is critical for spectral index imag- 
ing (iLobanovIl 99611 1998bl) and synchrotron tu rnover frequency 
imagmg ("Lobano v et al J 1 99 7;'Loban ovl 1 998al) . The effect of the 
core shift is illustrated by Fig. |7}{i]for radio spectral index im- 
ages of the jet in 3C 120. The spectral index image (Fig.|7]) and 
the spectral index profile along the jet ridge line (Fig.|8]l obtained 
without accounting for the core shift can clearly result in mis- 
leading conclusions about the physical conditions in the jet flow. 
Therefore, the alignment of VLBI images should be given proper 
consideration in all studies involving multi-frequency measure- 
ments. 

At present, the S / X band VLBI astrometry pr o grams 
(e.g.. lMaetal.1 1 1998; Fe vetal.1 120 04; Beasley et al.' 12002 



Fonialont etaill2003; Petrov et all 12005. ,2006; Kovalev et all 
2007h employ an implicit assumption that the source positions 
are the same for both bands (i.e., all objects are assumed to be 
point-like, optically thin radio sources). However, both the core 
shift and the source structure should be taken into account in 
order to achieve a higher positional accuracy of VLBI astrom- 
etry. Th i s issue regarding the source structure was discussed by 
I Chariot! (1 1 9901) . The likely time variability of the core shift (due 
to flaring activity of the nucleus. Sect. 13.21 ) may lead to even 
larger errors in positional accuracy of repeated astrometric VLBI 
experiments (e.g., the RDV project) and must be also taken into 
account. 



3.4. Radio-optical alignment of astrometric positions 

The core shift issue becomes an even more important fac- 
tor when the radio reference frame needs to be connected to 
an optical reference frame. So far the link is based on the 
study of optically bright radio-emitting stars which are seen 
both by Hipparco s and VLBI (see, e.g., Lestrade et al. 1995; 
iJohnston & de Veg t 1999; Ros 2005). Future accurate alignment 
of the frames has to rely on using compact radio sour ces in dis- 
tant quasars (e.g.. lFev et al.ll200lHBourda et al.ll2007l) . 

Below, we discuss the alignment problem for compact ex- 
tragalactic radio sources. We assume that the dominating com- 
ponent in both the radio and optical bands is the synchrotron 
self-absorbed compact jet origin (core). Broad-band modeling of 
blazar spectral energ y distribu tion supports this hypothesis (see, 
e.g., recent review by iB ottcheijl2007h . High-resolution VLBI ob- 
servations of nearby AGN imply that the jet is formed and emit- 
ting in the radio akeady at dist ances of < 1000 gravitational 
radii from the central engine (e.g. lJunor et al.l[T999l : lKadler et alJ 
[2004). Thus the physical offset between the jet base and the cen- 
tral nucleus is much smaller than the positional shift due to opac- 
ity in the jet (the latter can be larger than 1 pc). This implies 



o 




4 6 8 10 
Distance (mas) 

Fig. 8. Spectral index profiles along the jet ridge line obtained 
from the three different image alignments illustrated by the spec- 
tral index images of 3C 120 presented in Fig.|7] The dotted line 
corresponds to the alignment at the phase centers. The dashed 
line corresponds to the alignment at the model-fitted core posi- 
tions. The solid line corresponds to the alignment at the model- 
fitted positions of optically thin jet components. 



that the offset between radio and optical positions of reference 
quasars will be dominated by the core shift even if the optical 
emission comes from the accretion disk around the central nu- 
cleus. 

The magnitude of the core shift, Ar, between two arbi- 
trary frequencies vi and V2 (vi > V2) caused by synchrotron 
self-absorption can be predicted for an obje ct with known sy n- 
chrotron luminosity, Lsyn, of the compact jet (lLobanovll998bl) : 



Ar 
[mas] 



4.56 ■ 10" 



2/3 



~([pc]) ( 



1 +z 



r20O Wl -ySj COS0j)0 
i-syn \ / Vi V2/(Vi - V2) ] 



[erg/s] 



[GHz] 



(3) 



where Dl is the luminosity distance to the object, /3j and Fj are 
the jet speed and Lorentz factor, flj is the angle between the jet 
velocity vector and the line of sight, 0o is the jet opening angle 
and = ln(rn,ax/''min) describes the extent, (r,nin, ''max), of the 
compact jet. 

If not corrected for, the core shift will introduce an additional 
additive error factor in the alignment of the radio and optical 
reference frames. The magnitude of this error can be estimated 
from the following considerations. If the two reference frames 
are aligned using a sample of AGN with compact radio jets, the 
expected mean offset of the source position caused by the core 
shift can be calculated from equation ([3j. For this calculation, the 
sample averages (Dl), (Lj) and (Lsyn) can be applied. A conser- 
vative estimate of the jet opening angle 0o ^ 1 /L^ can be made, 
assuming that the jet has a Mach nu mber M\ ~ Fj. A typical 
compact jet extent of rn,ax/'"min = 100 (iLobanovll 1 998bl) can be 
adopted for all objects. It should also be noted that this parameter 
contributes to the result of the calculation only logarithmically. 

If the sources are selected from a flux density limited sam- 
ple, the distributio n of their orientations will be affected by the 
Doppler boosting jVermeulen & Cohenlll994t iKellermann et al.l 
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Table 2. Theoretical predictions of expected core shifts values 





Slim 

[deg] 


^med 

[deg] 


8.6 GHz V 
2.3 GHz' 

[mas] 


/A ^6000 A \ 
8.6 GHz' 

[mas] 


[10-*''erg/s] 


5 


107.4 


37.4 


0.13 


0.05 


6.2 


10 


69.2 


23.2 


0.19 


0.07 


3.5 


15 


55.2 


17.8 


0.23 


0.08 


2.6 


20 


47.4 


14.8 


0.26 


0.10 


2.2 


25 


42.1 


12.9 


0.28 


0.11 


2.0 


30 


38.3 


11.5 


0.30 


0.12 


1.8 



Column designation: Fj 

ing angle for detection; 



— jet Lorentz factor; t'li,,, — maximum view- 
'fmcd — median viewing angle for a sam- 



ple with random orientation; A ''2 jghz — average shift expected be- 
tween 2.3 GHz and 8.6 GHz; Arf^^_^ — average shift expected be- 
tween 8.6 GHz and 6000 A. AH parameters are calculated for a sample 
with (Lsyn) = lO-^^erg/s, (Dl) = 3.32 Gpc (z = 1). The last column 
lists the mean synchrotron luminosities (L'^^} required to reconcile the 
predicted (AtjI™^) with the median value of the core shift (0.44 mas) 
measured in the data presented in this paper. 



l2004tlCohen et alJl2007h . with the resulting maximum viewing 
angle given by: 



arccos 



(rf - 1)1/2 



(4) 



where To is the intrinsic brightness temperature of the jet, Tn^ 
is the minimum detectable brightness temperature in the obser- 
vations, s - 1/(3 - a) and g is t he spectral index of the ra- 
dio emission (iLobanov et al.l 2OO0 ) . In the cores of the VLBI 
iets. To ~ 10" K (ILobanov et alJ 120001: iKovalev et all l2005t 



iHoman et al.l |2006 |) , and typical VLB! snapshot observations 
reach Ty,^ ~ lO** K ('LobanovetiDllOOO'). 

With the assumptions described above, equation (O can be 
used to calculate the average core shift expected from observa- 
tions of a flux density limited sample of extragalactic objects: 



(Ar) 
[mas] 



1.65-10-21 /(Dm>\ 7 (isyn) 



Vl V2/(Vl 



[pc] 



[erg/s] 



2/3 



[GHz] 



(5) 



where Dm - DhIO- + z) is the metric (or "proper motion") 
distance to the object and evaluates the term [sin6'j/j6j(l - 
ySj cos ^j)]^/^ at a "median" angle of the sample, ^med, determined 
from the condition ^gA9 = 2 j^""^ dO. 

Table |2] gives average core position shifts expected between 
2.3 GHz and 8.6 GHz and between 8.6 GHz a nd 6000 A for an 
average synchrotron luminosity of erg s"' (lLobanovll998bl) 
and a quasar sample peaking at a redshift z = 1 . The shifts are 
calculated for several different values of (Fj) ranging from 5 to 
30. The average 2.3-8.4 GHz core shift of 0.3 mas predicted for 
(Fj) = 30 is close to the mean value of 0.44 mas measured from 
our data. The direct comparison of these two values is not strictly 
valid, as our measurements are not done for a complete sample. 
We can note however that this discrepancy can be easily rec- 
onciled by allowing for small variations of the model parame- 
ters (i.e. varying the mean luminosity of the sample, as illus- 
trated in Table|2]). With this assumption, the maximum measured 
value of the 2.3-8.4GHz core shift can be readily explained by a 
one sigma deviation from the mean in a sample of compact jets 



with a Gaussian distribution of Lorentz factors characterized by 
<Fj) = 30 and err = 10. 

It is immediately seen from Table |2] that the average shifts 
between the radio and optical bands are comparable to the po- 
sitional accuracy of VLBI and GAIA, and significantly exceed 
that of SIM. These numbers, albeit not definite and assumption- 
dependent, imply that the core shifts should be carefully investi- 
gated, and corrected for, before attempting to make a connection 
between the ICRF and the reference frames that will be produced 
by GAIA and SIM. 

To achieve this goal, a coordinated program should be estab- 
lished, aiming at defining a Primary Reference Sample (PRS) 
of extragalactic objects to be used for radio-optical reference 
frame alignment and investigating the magnitude and variabil- 
ity of the core shift in this sample. This would require making 
multi-frequency and multi-epoch VLBI observations of the PRS 
in order to measure the core shifts and monitor their variabil- 
ity for a period of time comparable with the mission duration 
of GAIA and SIM. Additional VLBI observations of the sample 
would have to be done contemporaneously with the optical mea- 
surements. 

The accuracy of the determination of the optical position of 
the core from radio measurements can be assessed in the follow- 
ing way. Let there be - 1 measurements of core shift between 
radio frequencies (vi , v^). Then the core position at an arbi- 
trary frequency can be determined from the core offset measure, 
Q.„ (see definition in Lobanov 1998b), provided that Q.„ ^ const 
for all measurements (this ensures that the accuracy of extrapola- 
tion does not depend on the actual wavelength to which the core 
position is extrapolated). The error, cry, can be estimated from 



A^- 1 



A'- 



N I .7 



1/2 



(6) 



where cri is the positional measurement accuracy at the lowest 
frequency, vi, in the data set. For typical VLBI measurements 
(SNRcore ~ 500, SNRjet ~ 100) spanning four frequencies from 
8 GHz to 43 GHz, we estimate cr^ k 0.05 mas, with referenc- 
ing the core position to an optically thin feature in the jet (mea- 
surements at 2 GHz and 5 GHz may still be required for moni- 
toring Q,v at lower frequencies and verifying that it is not vari- 
able). For relative astrometry VLBI measurements with phase 
referencing, we obtain a-y ^ 0.04 mas, assuming that measure- 
ments are done at four frequencies between 5 GHz and 22 GHz 
and respective positional errors are the typical ones reported by 
Fomalont (2005) for the case of one phase reference calibrator 
at a distance of 3° from the target. These estimates indicate that 
extrapolating core shifts to optical wavelengths can be done with 
sufficient accuracy for astrometric applications. 

If the core shift effect is not corrected for, reaching the accu- 
racy goal for the radio-optical alignment would require increas- 
ing the number of reference sources by some factor Using the 
average radio-optical core shift in Table|2] we estimate this factor 
to be about 3.5 and 4.5 for the alignment of the ICRF reference 
frame with the reference frames of GAIA and SIM, respectively. 
These values may be reduced by about 20 percent, if individual 
radio-optical position offsets are wei ghted appropriately w ith re- 
spect to the radio jet direction (e.g.. Kova lev et ani2005h along 
which the core is expected to shift. But even then it would re- 
quire a substantial increase in the number of reference objects 
required for reaching the sufficient accuracy of the alignment. 
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4. Summary 

Measurements of the frequency-dependent shift of the parsec- 
scale jet cores in AGN are reported for 29 bright extragalactic 
radio sources. It is shown that the shift can be as high as 1.4 mas 
between 2.3 and 8.6 GHz. If not taken into account, such core 
shifts could influence and corrupt both astrophysical (e.g., spec- 
tral imaging and Faraday rotation imaging) and astrometric stud- 
ies. It is shown that nuclear flares can result in temporal variabil- 
ity of the core shift. Analysis of this effect can shed light on 
the physics of the flares. The core shifts are likely to influence 
the positional accuracy of the radio reference frame and pose 
problems for connecting radio and optical reference frames. The 
next steps in this study are to measure core shifts in a complete 
flux-density-limited sample of bright extragalactic radio sources, 
to estimate physical parameters from the measured shifts and 
to look for their temporal variations. We plan to achieve these 
goals by combining together further analysis of RDV VLBA 
S/X datasets, four-frecmency (8. 1, 8.4, 12, & 15 GHz) MOJAVE- 
2 VLBA observation^ in 2006, and dedicated multi-frequency 
1.4-15 GHz VLBA observations in 2007 of more than half of 
the objects in Table [1] 

We have estimated from theory an average shift between the 
radio (4 cm) and optical (6000 A) bands to be approximately 
0.1 mas for a complete sample of radio selected AGN. The ro- 
bustness of this prediction is supported by the fact that the same 
method gives an average core shift of 0.2 to 0.3 mas between 
2.3 and 8.6 GHz, which agrees well with the median shift of 
0.44 mas reported by us here for a non-complete sample of 29 
AGN. 

The estimated radio-optical core shift is comparable to the 
positional accuracy of GAIA and significantly exceeds that of 
SIM. It implies that the core shift effect should be carefully in- 
vestigated, and corrected for, in order to align accurately the ra- 
dio and optical positions. We suggest two possible approaches 
to tie the radio and optical reference frames together. 1) In the 
first approach, multi-frequency VLBI measurements can be used 
for calculating the projected optical positions, assuming that the 
radio and optical emission regions are both dominated by a spa- 
tially compact component marginally resolved with VLBI and 
SIM and point-like for GAIA. The core shift measurements of 
the objects are required to be made quasi-simultaneously with 
astrometric measurements of their positions. The discrepancies 
between the measured optical and radio positions can then be 
corrected for the predicted shifts, and the subsequent alignment 
of the radio and optical reference frames can be done using stan- 
dard procedures. 2) A more conservative approach may also be 
applied, by employing the VLBI observations to identify, and 
include in the Primary Reference Sample, only those quasars in 
which no significant core shift has been detected in multi-epoch 
experiments. Either of the two approaches should lead to sub- 
stantial improvements of the accuracy of the radio-optical posi- 
tion alignment. 
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